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1 Introduction 

Holographic superconductors have recently received a new wave of attention. It orig¬ 
inated from several attempts [1-7] to provide a holographic description of systems 
which resemble more of the real-world superconductors. One of the essential fea¬ 
tures of the original holographic superconductor proposal of [8, 9] is that it describes 
the system which exists in two states: a superconducting state which has a non¬ 
vanishing charge condensate, and a normal state which is a perfect conductor. As a 
direct conseguence, already in the normal phase the static electric response, namely 
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the DC conductivity {u = 0), is infinite. This is a straightforward consequence of 
the translational invariance of the boundary field theory, which leads to the fact that 
the charge carriers do not dissipate their momentum, and accelerate freely under an 
applied external electric field. Therefore one is motivated to introduce momentum 
dissipation into the holographic framework, breaking the translational invariance of 
the dual field theory. It is definitely interesting to construct a holographic supercon¬ 
ductor on top of such dissipative backgrounds which is indeed going to have a finite 
DC conductivity in the normal phase, clearly distinguishable from the infinite one in 
the superconducting phase. 

One efficient method to implement such a feature relies on the possibility of 
breaking diffeomorphism invariance in the bulk via giving the graviton a mass, as it 
has been proposed in [10]. It is very convenient to recast these Lorentz symmetry 
violating massive gravity theories into a covariant form introducing the Stueckelberg 
fields, namely the extra degrees of freedom appearing as a consequence of breaking 
of the diffeomorphism symmetry (see [11] for more details). 

In the context of applied holography this construction was analyzed for the first 
time in [12] where momentum dissipation in the field theory was achieved by switch¬ 
ing on neutral scalar operators depending linearly on the spatial coordinates of the 
boundary. These scalar fields on the boundary source the neutral scalar fields in 
the bulk. The resulting bulk system describes a holographic dual of the field theory 
with broken translational symmetry. Such a system possesses a finite DC conductiv¬ 
ity [12]h 

The original idea of [12] has been put in a broader context in [17], where the most 
general form for the Lagrangian of the neutral scalars has been introduced^. This 
Lagrangian is weakly constrained by the consistency conditions in the bulk, which 
avoid ghost excitations and gradient instabilities [17]. It turns out that imposing 
physical consistency of the theory still leaves enough freedom to construct models, 
which exhibit new non-trivial features. 

To be more specific, one can build models which possess the following attractive 
properties. The first one is an increase of conductivity as a function of temperature, 
for temperatures lower than a certain critical value Tq, 

r<r„, (1,1) 

^See [13-16] for further studies about Massive Gravity as an effective description for Momentum 
Dissipation. 

more restrictive generalization has been analyzed in [18]. 
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This property bears a resemblance to an insulating behavior, with the population of 
the conducting energy band depleting upon lowering the temperature. Still, it awaits 
a better understanding, because of an essentially non-vanishing value of the DC con¬ 
ductivity at zero temperature. We refer to the state (1.1) as pseudo-insulating. The 
second new feature of the model is an appearance of an extra structure in the optical 
conductivity. For temperatures lower than a certain critical value T', there appears 
a peak in the optical conductivity, signaling a new long-lived collective propagating 
excitation of the charge carriers^. 

This paper is based on the idea to generalize the construction of [5, 6] to the 
more generic effective models for momentum dissipative systems, proposed in [17]. 
The main questions which we aim to answer are the following: 

1. Can one construct a model of holographic superconductor which is separated 
by the lines of the second order phase transition from the normal metallic phase and 
the normal pseudo-insulating phase (1.1)? 

2. Does the peak in the optical conductivity of [17] continue to exist in the 
superconducting phase^? 

We have found that the answers are: 

1. Yes, by combining the idea of [17] with the setting of a holographic super¬ 
conductor one can obtain a system with a rich phase diagram where three different 
phases are present: superconductor, metal, and pseudo-insulator. 

2. The peak in optical conductivity continues to exist in the superconducting 
phase, as the temperature is lowered below a critical temperature Tc of the super¬ 
conducting phase transition. However, at a certain temperature T = T" the peak 
disappears. 

Furthermore, we attempt to construct a holographic model with the phase di¬ 
agram containing a superconducting state inside a dome-shaped region, (see, e.g., 
[21] for the discussion about dome phase diagrams in condensed matter theory and 
figure 1 for a sketch of the realistic situation®). The most successful result would 

^It is really tempting to make a comparison to polaron physics, see, e.g., [19]. 

^See also [20], where non-trivial structure has been observed in the optical conductivity of a 
holographic superconductor. 

®Note that in nature the axes are temperature and doping while in our case doping is replaced 
by disorder-strength; in this sense we do not aim to construct an holographic dual of a real dome¬ 
shaped phase diagram typical of high-Tc superconductors. 
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Figure 1. Schematic phase diagram of a real cuprate high-Tc superconductor. 

be to have a superconducting dome, separated from an insulating normal state at 
smaller values of the disorder strength parameter, and a metal normal phase at its 
larger values. In hgure 1 we provide a schematic sketch of what we would like to 
approach, the phase diagram for High-Tc superconductors. 

We will demonstrate that implementing the momentum dissipation models of [17] 
in the holographic superconductor framework can indeed lead to a superconducting 
dome, located between pseudo-insulating and metallic phases. However, it appears 
that such models are too restricted to describe superconducting dome with realistic 
critical temperature of the superconducting phase transition. We have found that 
the critical temperature of the dome Tc{a), where a is the magnitude of the trans¬ 
lational symmetry breaking, is bounded from above by a small number (in units of 
charge density), of the order of 10“®. This makes the numerical calculation at hnite 
temperature hopeless. 

Nevertheless, at zero temperature it is possible to have analytical control of the 
SC instability through the BF bound reasonings and show the existence of a super¬ 
conducting dome. In the Discussion section 8 we provide a few ideas to generalize 
our model, which might be useful to obtain a superconducting dome with a reason¬ 
ably higher values of the critical temperature. We will be considering charged black 
brane backgrounds with the neutral scalar helds having vacuum prohles, depending 
linearly on the spatial coordinates: 


(t)^ = ax, (f)^ = ay. ( 1 . 2 ) 

This conhguration (1.2) breaks translational symmetry (and Lorentz invariance) of 
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the boundary field theory but keeps untouched energy conservation. Within this 
choice we are going to retain homogeneity and rotational invariance. It would be 
interesting to reproduce the same sort of computations in an anisotropic setup as in 
[7]. Besides the parameter a, describing the magnitude of the translational symme¬ 
try breaking, we will also introduce another parameter m, which will be primarily 
important in the models with non-linear action for the neutral scalars (1.2). 

We consider the system at a finite charge density, which corresponds holograph¬ 
ically to the time-like component of the U{1) gauge field having a non-trivial radial 
profile in the bulk, At{u). The charged scalar -0 is dual to the condensate O of charge 
carriers. When the v.e.v. of the condensate is non-vanishing, (O) ^ 0, the system is 
in a superconducting phase. This corresponds holographically to a non-trivial config¬ 
uration 'ip{u) in the bulk, with the vanishing source coefficient of the near-boundary 
expansion of the 'ipiu) [8, 9]. 

We will study various superconducting systems, distinguished by the choice of 
the Lagrangian V{X) for the neutral scalar fields, where 

v = (1.3) 

and L is the radius of AdS. In this paper we will be mostly interested in the following 
models®: 


model 1 : 

1-(A) = — , 

2 

(1.4) 

model 2 ; 

V{X) = A + X® 

(1.6) 

model 3 : 

wAf 

1-(X) = —, X^l 

2 

(1.6) 

The model (1.4) gives the simplest 

way to describe the fields 0^ 

and has been pro- 


posed in [12]. 

We will argue that already in the simple case of (1.4) it is possible to have a 
superconducting dome. We will demonstrate this analytically at zero temperature. 
Interestingly, the dome is achieved for the scaling dimension A and the charge q of 
the scalar -0, restricted to the small vicinity of the “dome” point, which we have 
found to be 

(Ad, Qd) = (2.74, 0.6). (1.7) 

®In models 1.4 and 1.6 we introduced the prefactor of 1/m? into the definition of V{X). Such 
change of notation will render a to be the only translational symmetry breaking parameter in the 
models 1.4 and 1.6. In model 1.5 instead both m and a are independent parameters and we decided 
to avoid any rescaling in the definition oiV{X). 
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In this case the superconducting dome exists in the middle of a normal metallic phase 
(the model (1.4) does not allow an insulating phase). 

We will show that the model with the non-linear Lagrangian (1.5) also possesses 
the superconducting dome near the point (1.7). In this case it is possible to engineer 
a model where the dome is separated from metallic phase at larger values of the 
translational symmetry breaking parameter m, and from a pseudo-insulating phase 
at smaller values of m. This situation is qualitatively the closest one to the actual 
real phase diagram for High-Tc superconductors. It is important no notice that our 
dome is constructed dialing the disorder-strentgh parameter of the theory, while the 
actual dome in High-Tc Superconductors depends on the doping' of the material. 
We are not aware of experimental phase diagrams where the SC dome occurs as 
function of increasing disorder-strength. 

To support our statement about the superconducting dome with such a small 
critical temperature T^, we will calculate numerically the dependence of the critical 
temperature for the models (1.4), (1.5), on the scaling dimension A and the charge 
q. We will show that as the (A, q) approach the dome point (1.7), the critical tem¬ 
perature quickly declines. 

The rest of this paper is organized as follows. In the next Section 2 we set up the 
model which we will be studying in this paper. We consider the general Lagrangian 
V{X) for the massless neutral scalar fields. In Section 3 we review the properties of 
the normal phase solution. In Section 4 we study the conditions for its instability 
towards formation of a non-trivial profile of scalar hair. From the field theory point 
of view this corresponds to a superconducting phase transition. In Section 5 we 
focus on the features of the broken phase, the condensate and the grand potential, 
demonstrating explicitly the second order phase transition at T = Tc. In Section 6 
we study the optical conductivity in the normal and superconducting phases. In 
Section 7 we describe the way to construct a superconducting dome in the middle 
of a metallic phase, for the model (1.4), and between pseudo-insulating and metallic 
phases, for the (1.5). We discuss our results in Section 8. Appendix A contains 
further details about the calculations of the condensate and the grand potential. 
Appendix B is dedicated to derivation of the on-shell action for bulk fluctuations, 
which are holographically dual to current and momentum operators on the boundary. 

Tf our parameter was the doping of the material it would affect the charge density of the 
boundary theory and this does not happen in our model. 
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2 Setting up the model 

In this section we introduce the model, which we will be studying in this paper. 
We begin by writing down the action and equations of motion of the bulk theory. 
We proceed by deriving equations of motion for the general ansatz, describing the 
charged black brane geometry, with linearly-dependent sources for 0 ^, and radially 
dependent charged scalar Then we will review the normal-state solution of the 

model, which has a trivial charged scalar field profile = 0 . 

2.1 Action and equations of motion 

The total action of our model is : 


I — h + h + h ^ 


( 2 . 1 ) 


where we have denoted the Einstein-Maxwell terms A, the neutral scalar terms / 2 , 
and the charged scalar terms J 3 ; 


R-2A- — 


h = j 

R = -2m^ j {X) , 

h = - [ d<^+^x^g + + |^| 2 ) 


( 2 . 2 ) 


We have inserted an additional coupling m? in front of the potential V{X) which is 
going to be redundant for the monomial cases 1.4 and 1.6 where we decided in fact 
to reabsorb it into the definition of V {X). In this way for those cases we are left with 
just one parameter a which is going to represent the disorder-strength in the system. 
In the case of the polinomial potential 1.5 is going to be an independent parameter 
in addition to a. We have introduced an extra coupling k, between the charged scalar 
-0 and the neutral scalars (j)A In this paper we will be mostly considering k = 0, and 
comment on the models with non-vanishing k in the discussion section 8 . We have 
defined 

X = (2.3) 

We denote = (d^ — iq A^)'ip to be the standard covariant derivative of the scalar 
'ijj with the charge q. We fix the cosmological constant to be A = —3/L^. In this 
paper we will consider 4-dimensional bulk, d = 3. 
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The equations of motion following from the action I read® : 


- \D^P\^ - (M^ + nHm^ - 2mV^ 

= (^m‘^V + d^4>'d.ct)^ + yF^A^"/ + \{D^^pD,r + D^i^D^r) 

^ d, F'^^) - - ^D’^r) = 0 (2.4) 


1 


D^[^g D^i^) - {M^ + KH)i: = 


df, (y^ (2m^ V + kH\'iI:\‘^) = 


where the dot stands for a derivative w.r.t. 


0 

0 , 


(2.5) 

2.2 Background 

We consider the following black brane ansatz for the background: 

ds^ = + \{dx^ + dy^) + ], 

\ u^j[u) j 

(I)^ = a5^x\ I,i = x,y. ( 2 . 6 ) 

A = At{u)du , V’ = ■ 


The 0^ scalars have prohles linear in the spatial coordinates x, y oi the boundary. 
They effectively describe momentum dissipation mechanisms in the boundary held 
theory, making the DC conductivity of the theory hnite® [12]. We will take 0 to 
be real-valued, since due to the u component of Maxwell equations the phase of the 
complex held 0 is a constant. We are looking for charged black brane solutions with 
a scalar hair where Uh is the position of the horizon, and the boundary is located at 
M = 0. We allow for non-trivial x{'^) because we want to have in general a non-trivial 
'0 (m). If 0 = 0, then x = 0- 

®When V{X) = Xf2m‘^, we recover the equations of [6]. 

®These fields are dual to marginal scalar operators whose sources explicitly break translational 
symmetry. Exploiting the shift invariance for these operators it is possible to retain the homogeinity 
of the background such that the metric and the charged scalar (and as a conseguence the stress 
tensor and the SC order parameter) do not depend on the spatial coordinates at all. 






(2.7) 


The resulting equations of motion read: 
u Af -0^ 




p + u'tp'‘^ = 0 


P 


1,2 2 /' e^u^Af nL^H'ip^ e^q'^A'fp" 

2m^L‘^V 2AL2 6 
2q‘^Ati>^ X' 

~PPr~ 2 


- PA[ - A'l = 0 


;// / 2 f xP ,, [e’^q^Al kH ^ 


= 0 


The Hawking temperature of the black brane (2.6) is given by: 

rr, f'{Uh) 

1 = -^- e 2 . 

dvr 

Using eqs. (2.7)-(2.10), the temperature can be written as: 

_ X 

T = - (-12 + WlV + 2(M2 + KH)L‘^'tp‘^ + e^ulA'P . 

IQnUh 

with all the helds evaluated at the horizon Uh- 


( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


2.3 Normal phase 

In the case of a non-trivial condensate 'ip{u) it is in general impossible to solve the 
background equations of motion (2.7)-(2.10) analytically. However, when ^lJ{u) = 0, 
the solution is known [17]. 

From now on we will £x the coupling k to zero, 


K = 0 . 

The resulting normal phase background is given by: 


(2.13) 


^(7 = 0, x(m) = 0, 


At{u)=^-up, p = puh, 


PP 



( 3 m^L‘^V{aP^) 

V ^ 



dy 


(2.14) 

(2.15) 

(2.16) 


Due to (2.12) the temperature in the normal state reads: 

T = - (-12 + dm^LV + uy) . (2.17) 

IbTTUji 

All the features of this normal phase solution are going to be reviewed in detail in 
the following section. 


9 



















UJ 


Figure 2. Left: Optical conductivity in the normal phase for (1.5) (a = ^/2, m? = 0.025, 


p = 1) with temperature running from T = 0.04 (black line) to T = 1 (red line). Right: 
DC conductivity as a function of temperature for the same model and parameters. 

3 Normal phase features 

As suggested in [17], for models with a specific choice of the Lagrangian V{X), the 
solution exhibits various interesting properties. 

Using the membrane paradigm the DC part {u = 0) of the optical conductivity can 
be computed analytically [15] and for a generic Lagrangian V{X) it is given by [17]: 



(3.1) 


The DC conductivity consists of two parts: 


(3.2) 


^DC ^pair T ^dissipation ) 


which is a generic holographic feature The first one Cpair is due to pair creation in 
the background, and it is present even at zero charge density [22]. It corresponds 
exactly to the probe limit result. It is temperature independent, and therefore is 
always present (unless we introduce a dilaton held) as an offset in the value of aoc, 
leading to aociT = 0) 7 ^ 0. The second term a dissipation is really the one dealing with 
dissipative mechanism, and it can be thought as the strongly coupled analogue of 
the Drude formula for the conductivity. In the limit of zero translational symmetry 
breaking parameter m, this second term gives rise to the inhnite DC conductivity, 
typical for backgrounds preserving translational symmetry, such as the AdS Reissner- 
Nordstrom black brane case. Due to the freedom of choice of the Lagrangian V{X) 
this solution can be either a metal or a pseudo-insulator and can provide a transition 
between the two phases (see hgure 2). The pseudo-insulator phase is characterized 
by the conductivity, declining at smaller temperatures, da/dT > 0, for T < Tq, but 
reaching a non-vanishing value at T = 0 (which is the reason why we are not calling 
it an insulating phase)The transition between the two phases is provided by the 

One easy way to enable ctdcIF = 0) = 0 is adding a dilaton field to the action [23], which 
allows to get a “real” insulating state. See also [24] for an alternative approach. 
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existence of a maximum in the DC conductivity as a function of temperature (see 
figure 2), at T = Tq, which gives a clear separation between two different regimes: 


da 

df 

da 

df 


< 0 , 

> 0 , 


T > To —)■ metal 
T < To —)■ pseudo-insulator 


(3.3) 

(3.4) 


The temperature Tq at which the metal-insulator transition happens can be obtained 
analytically, solving the following equation: 


da DC 

dufi 


YV{Y) = V{Y), Y = uW. 


(3.5) 


The metal-insulator transition in the behavior of the DC conductivity is related 
to a non-trivial structure in the optical conductivity, namely a weight transfer from a 
Drude peak into a localized new peak in the mid-frequency regime (see figure 2). This 
feature corresponds to an emerging collective propagating excitation of the charge 
carriers, whose nature is not completely clear yet. The phase diagram of this nor¬ 
mal phase is already rich and can give insights towards the interpretation about the 
various ingredients introduced into the model. In the case of the linear Lagrangian, 
which goes back to the original model [5], the parameters m and a are combined into 
ma, which can be interpreted as the strength of translational symmetry breaking. 
From the dual field theory point of view this is thought to be related to some sort of 
homogeneously distributed density of impurities, representing the disorder-strength 
in the material. 


In the case of a more general V{X), the m parameter keeps this kind of inter¬ 
pretation while the a one represents the strength of interactions of the neutral scalar 
sector. This reasoning is conhrmed by the study of the phase diagrams of the system 
(figure 3) which makes evident the difference between the two parameters. Indeed, 
while the m parameter, which we are going to interpret as the disorder-strength of 
our High-Tc superconductor, enhances the metallic phase, the a one clearly reduces 
the mobility of the electronic sector driving the system towards the pseudo-insulating 
phase. 


4 Superconducting instability 

In this section we will describe the instability conditions for the normal phase towards 
the development of a non-trivial prohle of the charged scalar held. This allows one 
to determine a line of the second order superconducting phase transition, Tc(a) (or 
Tc(m) for the model (1.5)), in the boundary held theory, with broken translational 
symmetry. 
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■ Pseudo-Insulator □ Metal 


■ Pseudo-Insulator □ Metal 



m'^2 a 

Figure 3. Region plots for the model (1.5) in the normal phase. We choose units where 
the density is p = 1. Here we have fixed a = 1 (left plot) and m = 1 (right plot). The blue 
region is pseudo-insulating, daoc/dT > 0, the green region is metallic, daoc/dT < 0. 


We start by considering the system at zero temperature, which we are able 
to study analytically. Then we proceed to studying the normal phase at a hnite 
temperature. Upon lowering the temperature, at a certain critical value T = Tc, 
the normal phase becomes unstable. This is the point of a superconducting phase 
transition. We construct numerically Tc as a function of the parameters A, q, a (or 
m), for the models with various V{X). 


4.1 Zero-temperature instability 


In the case of T = 0 the normal phase geometry interpolates between the AdS 4 ^ 
in the ultra-violet and the AdS 2 x in the infra-red. We can apply the known 
analytical calculation to study the stability of the normal phase towards formation 
of a non-trivial prohle of the scalar -0 [25]. 

Due to eq. (2.10), the effective mass Me// of the scalar -0 is given by: 


+ K H A^^ . (4.1) 


Notice that at the boundary the mass of the scalar is just M^ but at the horizon it 
gets an additional contribution. This is because near the horizon we have: 






L2 f"{uh){Uh - u) 


(4.2) 


at zero temperature. Due to (2.16), and the zero temperature T = 0 condition, with 
the temperature given by (2.17), we obtain: 


/"K) 


2 ^6 -|- L’^rn? 2U(w^a^) + u\a^V (m^q;^) j j 


(4.3) 
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Figure 4. Region and contour plots for the model (1.4) with linear potential for the 
neutral scalars (dimensional quantities are measured in units of the charge density p). We 
choose a = 2. The region of A, g, satisfying the IR instability condition (4.6) is shaded 
in grey. The red dot is centered around {qd,Ad) = (0.6, 2.74). These tuned (g, A) confine 
superconducting phase of the model (1.4) into a dome region, as we discuss in Section 7. 
Notice the proximity of the red dot to the boundary of the IR instability region, resulting 
in Tc{qd, A^) being very small. 


The normal phase is unstable towards formation of the scalar hair, if Me// violates 
the BF stability bound in the AdS 2 , namely: 

(4.4) 

In (4.4) we have denoted the AdS 2 radius as L 2 ^^ : 

Combining (4.1), (4.3), (4.5), the IR instability condition (4.4) finally reads^^ : 

L><0, (4.6) 


where we have defined the function D as: 


^=4 


^ {kH + M^) f f a^u\V — 2V j + 6 j — q^p 


. 2 . 2^4 


L‘^m? ia^ufy {cAu’f) — 2V (a^uDj + 6 


(4.7) 


^^In the usual RN case we have f"{uh) = ^' 2 , and we find the usual = 4- in d = 3. 
^^This formula agrees with [6] in the case of V{X) = and /c = 0. 


13 











q 


Figure 5. Region and contour plots for the model (1.5) with non-linear potential for the 
neutral scalars with a = 0.25, m = 4 (dimensional quantities are measured in units of 
the charge density p). The region of A, g, satisfying the IR instability condition (4.6) is 
shaded. The blue dot on the plot has coordinates (g^, A^) = (0.6, 2.74). These tuned 
(g. A) confine superconducting phase of the model (1.5) into a dome region, as we discuss 
in Section 7. Notice the proximity of the blue dot to the boundary of the IR instability 
region, resulting in Tc{qd, A^) being very small. 

For the practical calculations we will solve the equation T = 0, see (2.17), for the 
value of Uh, giving the position of the horizon of the extremal black brane, 

- 12+ uIp‘^+ AL‘^m'^V{ula‘^) =0. (4.8) 

We will measure all the dimensional quantities in units of p; both for zero temperature 
and hnite-temperature instability analyses the p can be scaled out. 

In hgure 4 we plot the IR instability region on the (A, q) plane, for the model 1, 
(1.4), with a = 2, as well as a few contour lines of the constant critical temperature. 
In hgure 5 we plot the IR instability region and several Tc = const curves on the 
(A, g) plane, for the model 2, (1.5), with a = 0.25, m = 4. Analogous plot for 
ordinary holographic superconductor can be found in [25]. Plot in the case of the 
linear V{X) model hrst appeared in [6]. 

4.2 Finite-temperature instability 

Consider the system at large temperature in a normal phase, which exists in a super¬ 
conducting phase at low temperatures. Therefore as we decrease the temperature. 
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Figure 6. Critical temperature as a function of a for the model (1.6). Left: q = 1, A = 2. 
Right: g = 3, A = 2. 



Figure 7. Critical temperature as a function of A for the model (1.6). Left (a = l, q = 2). 
Right: (a = 1, q = 0.6). 


at certain critical value the superconducting phase transition occurs. If is 
non-vanishing, then for T < the system is in a superconducting phase, with a 
non-trivial scalar condensate 

Recall that near the boundary the scalar held with mass M: 

M = ^VA(A-3), (4.9) 

behaves asymptotically as: 

’ (4A0) 

where ipi is the leading term, identihed as the source in the standard quantization. 

To hnd the value of we can look for an instability of the normal phase towards 
formation of the scalar held prohle [25, 26]. Near the second order phase transition 
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Figure 8. Critical temperature as a function of a for the potential V{X) = Xj^rr? + 
f5X^/2rin? for different choices /3. All the curves have a runaway behavior at a —)> oo, and 


only the shape depends on the value of (5. 

point T = Tc the value of V’ is small, and therefore one can neglect its backreaction 
on the geometry. The SC instability can be detected by looking at the motion of 
the QNMs of 'ip in the complex plane. To be more specihc, it corresponds to a 
QNM going to the upper half of a complex plane. Exactly at critical temperature 
we have a static mode at the origin of the complex plane, a; = 0, and the source at 
the boundary vanishes, '01 = 0- In fh® next section we will solve numerically the 
equations (2.7)-(2.10) for the whole background, and conhrm this explicitly. 

The scalar held is described by eq. (2.10), which in the normal phase becomes: 



(4.11) 


where f{u) is given by (2.16). To determine the critical temperature we need to 
hnd the highest temperature, at which there exists a solution to eq. (4.11), satisfying 
the "01 = 0 condition. In this case for T < Tc the system is in a superconducting 
state, with a non-vanishing condensate 02 - 

We are interested in the phases of the models (1.4)-(1.6) on the temperature- 
disorder strength plane. In hgure 6 we plot Tc{a) for the model 1, (1.4), and the 
model 3, (1.6), with N = 1/2, 2,3, for different values of the charge q. It is clear 
that when the power N in the potential V is higher, the critical temperature for 
the SC phase transition is smaller. One interesting behavior, which still lacks an 
interpretation, is the non-monotonic behavior of Tc as a function of a, which was 
already observed in the original model [5, 6] and still persists in more generic setups. 

In hgure 7 we plot Tc{A) for q = 0.6 and a = 1 for the model 1, (1.4), and 
the model 3, (1.6), with N = 1/2, 2,3. The Tc(A) curves explicitly show that the 
critical temperature quickly declines as A approaches the border of the IR instability 
region. It is further underlined how higher powers/non-linearities in the potential 
lead to deeper supression of the critical temperature. 
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Figure 9. Condensate for A = 2, g = 1, au/j = 0.5, mL = 1 model with V = z + z^, and 
the corresponding Grand Potential for the two phases. 

We also plot Tc(a) for the generalized model 2, (1.5), in figure 8 for various 
amounts of non-linearity f3 , showing again the same behavior of suppression of 
the superconducting phase at larger (3. 

5 Broken phase and phase diagram 

In this section we study superconducting phase and construct the phase diagram on 
the (m, T) plane of the model (1.5). We will conhrm existence of the second order 
phase transition between normal and superconducting phases by solving four equa¬ 
tions (2.7)-(2.10) for the fully backreacted background. Knowing the near-boundary 
asymptotic behavior of this solution, one can determine the grand potential of the 
superconducting phase, and compare it with the grand potential of the normal phase 
to corroborate the phase transition at T = T^. 

Running the numerical procedure described in details in appendix A, we were 
able to construct the condensate as a function of temperature In 

hgure 9 we provide the plot of the condensate, for the model (1.5) with A = 2, g = 1, 
auh = 0.5 {a in units of entropy density), mL = 1. There we also plot the grand 
potential for the broken and normal phases which conhrms the superconducting 
transition at T = Tc- 

The holographic prescription for the calculation of the grand potential is: 


n = -T \ogZ = TSe, 


(5.1) 


where iSg is a Euclidean on-shell action of the bulk theory. 

After some computations showed in details in appendix A we obtain: 



(5.2) 
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Figure 10. Left: The T = 0 contour plots of the IR instability lines, for the model (1.5) 
with a = 0.7 and various choices of m. Right: Phase diagram for the model (1.5) at the 
point (7 = 4, A = 3 with a = 0.7 . We have chosen units p = 1. 

where we have also used the area-law expression for the entropy: 

s = lL 

The grand potential is dually given by: 

n = PV = -^SE, 

lOTT 

where V is a volume of spatial region. In conclusion we obtain (denoting p 
the expected thermodynamic relation: 

P = e-TS - pp. (5.5) 

where the energy density of the system is given by: 

e = -272^2 . (5.6) 

We now have enough information to construct the full phase diagram of the non¬ 
linear model (1.5). In dgure 10 we plot the phase diagram of the model (1.5) with 
A = 3, g = 4, a = 0.7 (in units of p = 1). We see that the superconducting region 
can be connected smoothly to both a metallic phase and a pseudo-insulating phase. 

6 Optical conductivity 

Our main aim in this section is to see whether the non-trivial structure in the optical 
conductivity (see dgure 2), pointed out for the model (1.5) in the normal phase [17], 
persists to exist in the superconducting phase. 


(5.3) 

(5.4) 
= pL=) 
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6.1 Fluctuation equations 

In order to compute the optical conductivity, we study the fluctuations on top of the 
charged black brane background with spatially-dependent neutral scalars, as follows: 

htxiu) 


^ 9 tx(t, u, y) 

r+oo 

dudk i^t + ikv 

J —oo 

{2ny 

tx, y) 

r+oo 

dudk 

/ —OO 

( 2 vr )2 

6Ax{t,u,y) = 

r+oo 

dudk 

J — OO 

(27r)2 




( 6 . 1 ) 


We consider homogeneous perturbations defined by fc = 0 , for which it is consistent 
to put all the fluctuations, besides (6.1), to zero. In this section we also put L = 1. 
The equations for the perturbations read^^ : 


2./,2X gxyl' 

CLx — 0 


// f f x'\ / f 2 -0 


OJ 


-C = 0 


e" 


2 /' x' 2ua^V{u^a^) 

u f 2 V{u‘^a‘^) 


e + 




P 




ie^ a u 


( 6 . 2 ) 

(6.3) 

htx = 0 (6.4) 


One can eliminate from the second equation (6.3) right away, and substitute 
it into equations for ax and ^ [12]. It is then convenient to perform the following 
redefinition: 

C(m) = . (6.5) 

tuj au^ 

and reduce the problem to a 2 x 2 system: 


e-^/O +e 
Pfe- 


2Pr^ eHu:^-PfA?y ^ a.+2m^aVe-tl/A;C=0 (6.6) 




f 


1 / 


e ^VC) 1 T P-^t^x + 


P V? 


2rrype ) C = 0 , 


/ 


(6.7) 


which in the normal phase agrees with the equations, derived in [17]. 

6.2 Superconducting phase 

In order to extract the optical conductivity of the system we first derive the on-shell 
action for the fluctuations. We leave the technical steps for appendix B while here 
we just quote the result: 


aP (cu) 


iL = J M i 


( 6 . 8 ) 
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In the case V{X) = these equations reduce to the fluctuation equations obtained in [6]. 
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where C^= Zju and we have dehned the matrix M. to be: 


M 


1 

0 

V 


\ 

0 

yi-2V2+^^^ ^ 




and expanded the fluctuations near the boundary u = 0 as: 


(u, oS) = (ca) + (ca) u , 

Z{u,u) = -4 {u,u ), 

luau 

Z{u,u}) = Z^^\uj) + Z^^\u})u. 


(6.9) 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 


We solve two coupled fluctuation equations (6.6), (6.7) numerically, for two inde¬ 
pendent sets of initial conditions which satisfy the infalling behavior near the horizon 
[27]^^. Due to linearity of the fluctuation equations (6.6), (6.7), the precise choice of 
the two sets of initial conditions is not important, and one can check that correlation 
matrix does not depend on it. For example, let us choose: 

{uh-u)7^^'^^ (6.13) 

= (_\) {uh - . (6.14) 


Near the boundary the fields behave as: 


ax^ 

Zb) 


Af 


+ \bU)]^^ ^' = 1 ^ 2 . 


(6.15) 


and we can assemble the matrices of leading and subleading coefficients: 




/lb)\ / r(1) d{2) 


.( 1 ) 4(2) 


B = 


0(1) d ( 2 ) 
Bz Bz 


(6.16) 


We collect the entries of the matrices (6.16) by integrating the equations for the fluc¬ 
tuations numerically and extracting the asymptotic behavior using (6.15). Knowing 
(6.16) and (6.9), we can Anally calculate the correlation matrix: 


g = MBA -'. 


(6.17) 
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See also [28] for an example of calculation of the correlation matrix in a different system of two 


coupled fluctuations. 
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Figure 11. The AC conductivity for the model (1.5) with a = y/2 (in units of l/u/j), 
m?L?‘ = 0.025, g = 4, and A = 2. Black line is at the temperature, slightly be¬ 
low the corresponding critical temperature ~ 0.16, and matches the result of 

the normal phase calculation at T = Tc- Red, blue, orange and green lines are for 
Tj = 0.15,0.12,0.09,0.06, respectively. Notice that as we decrease the temperature, 
between blue and orange line, the peak in the imaginary part of the AC condnctivity dis¬ 
appears. We call the corresponding critical temperature T”j ~ 0.1. We also provide 
the condensate as a function of temperature and mark the points where we calculated the 
AC conductivity. 

Finally from the correlation matrix (6.17), it is straightforward to hnd the AC 
conductivity in superconducting phase, using the Kubo formula: 

( 6 . 18 ) 

In hgure 11 we plot the AC conductivity for the model (1.5) with the non-linear 
Lagrangian for the neutral scalars, for A = 2, g = 4, a = y^, m?L?‘ = 0.025. We 
consider various values of temperature running from the normal phase to the super¬ 
conducting phase. The AC conductivity for the normal phase of the model (1.5) has 
hrst appeared in [17], where it has been shown that (temperatures are measured in 
units of square root of charge density) 
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1. For T > Tq (Tq ~ 0.46 (for the considered model) the system exhibits a 
metallic behavior, daoc/dT < 0. 

2. For T < Tf) the system exhibits an insnlating behavior, daoc/dT > 0. 

3. For T <T' < Tq, where T' ~ 0.35 (for the considered model) the non-trivial 
strnctnre in the AC condnctivity appears. To be more precise a Mid-freqnency peak 
shows up signaling a weight transfer mechanism and an emerging collective degree 
of freedom. 

These properties are illustrated in hgure 2. We checked that the sum rules for 
the optical conductivity are satished in both normal and broken phases. 

After we couple this model of [17], with the potential (1.5) for the neutral scalars, 
to the superconducting sector, more features appear. For the choice of parameters 
A = 2, q = A, auh = \/2, = 0.025 we continue to enumerate what happens as 

we decrease the temperature: 

4. At Tc ~ 0.16 (for the considered model) the second order phase transition 
occurs. The system lives in a superconducting phase, when T < Tc. 

5. At T = T" ~ 0.1 (for the considered model) the peak in the imaginary part 
of the AC conductivity disappears. The peak in the real part of the AC conductivity 
in superconducting phase gets smaller as the temperature is lowered and eventually 
disappears. 

These properties can be seen in figure 11. We will comment more on these fea¬ 
tures in discussion section 8. 

It would be very interesting to hnd the QNM excitations of the system in both 
normal and broken phase to have complete control on its transport properties and 
its collective excitations. We leave this topic for future studies. 

7 Dome of superconductivity 

In this section we describe how to construct a superconducting dome, by tuning 
the parameters of the model (1.5) with the non-linear Lagrangian for the neutral 
scalars. In nature, High-Tc superconductors exhibit a dome of superconductivity 
(see hgure 1) between insulating and metallic normal phases as a function of dialing 
the doping of the sample. We will construct a qualitatively similar behaviour but 
increasing the disorder-strength of the system. 

Due to limitations of our system we cannot construct an actual insulator, however the 


22 



q 



Figure 12. Searching for the dome for the model (1.5) with a = 0.25 at vanishing 
temperature. From the left graph we conclude that in order to have the superconducting 
dome we need to choose A and q from a small vicinity of the point (7.2). In the right 
graph we have verified explicitly the IR instability of the model with A = 2.75, q = 0.61, 
between two finite values of m. 


non-linear model (1.5) still allows to distinguish between two qualitatively different 
states of the normal phase, (3.3) and (3.4). 

The hrst observation is that when we decouple the translational-symmetry break¬ 
ing sector of the neutral scalar helds, by setting m = 0, we restore the framework of 
an ordinary holographic superconductor. Therefore, in order to conhne the super¬ 
conducting phase inside a dome, we need to make sure that the ordinary holographic 
superconductor exists in the normal state at any temperature. The way to achieve 
this is to make sure that the parameters A and q are such that the normal phase 
at T = 0 is stable. That is, we should have D > 0, where D is given by (4.7), with 
K = 0 and m = 0. The T = 0 IR stability condition therefore reduces to a well-known 
inequality, which reads: 

3 + 2A(A-3)-4g2 > 0. (7.1) 

Suppose now we stay on top of the T = 0 line on the (A, q) plane of ordinary 
holographic superconductor. The next step in engineering a model exhibiting a 
superconducting dome, is to restore a superconductor at a hnite value of m = mi, 
and then make sure that there is another value m 2 > mi, such that the system at 
m > m 2 is again in a normal phase. The procedure to search for the parameters 
which lead to the superconducting dome is the following. For the chosen value of a 
we plot D = 0 curves one the (A, q) plane, with the D given by (4.7), parametrized 
by various values of m. We search for the points (A, q) of intersection of two curves, 
corresponding to two different values of m. These values of m can be the boundaries 
of the dome region at T = 0. We then verify this explicitly by plotting the D for 
given a, A, q. In hgure 12 we plot the D = 0 curves for a = 0.25, on the (A, q) 
plane, and demonstrate explicitly that the dome requirement restricts us to consider 
a small sub-region on the (A, q) plane. In hgure 13 we repeat this for a = 0.5, and 
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Figure 13. Top: The D = 0 contours for the model (1.5) with a = 0.5; Bottom Left: 
BF bound violation for A = 2.728, q = 0.61; Bottom Right: Full Phase Diagram for the 
model. Green region is a normal pseudo-insulating phase, grey region is a normal metallic 
phase, red region is a superconducting phase. 

also plot the corresponding phase diagram. The superconducting phase is bounded 
from above by a small critical temperature, and is represented on the graph by a red 
interval. 

We have found that the requirement of having an interval of superconductivity 
[mi, m 2 ] at T = 0 is rather restrictive^®. We have found that on order to achieve the 
‘dome’ at a vanishing temperature we need to tune A and g to a small subregion of 
the region (7.1), centered around the point 

(Ad, qi) ~ (2.74, 0.6). (7.2) 

For such A and q we can engineer a model which, at T = 0, exists in a normal pseudo- 
insulating phase for m G [0, mi], in a superconducting phase for m G [mi, m 2 ], and 
a normal metallic phase for m > m 2 . 

Being more specific, it seems that there exists a lower bound for A below which no SC dome 
can be built within this model. It would be nice to understand better this bound. 
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Figure 14. Top: The D = 0 contours for the model (1.4); Bottom Left: BF bound 
violation for A = 2.745, q = 0.6; Bottom Right: Full Phase Diagram for the model. 
Orange region is a superconductor, lettuce region is a metal. 

The next step to construct the superconducting dome is to study the phase 
structure of the system at hnite temperature. To determine the boundary of the 
superconducting region, that is the line of the second oder superconducting phase 
transition, we can start in the normal phase, at larger values of temperature, and 
determine when it becomes unstable towards formation of the scalar hair. This 
procedure has been reviewed in Subsection 4.2. 

However, the point (7.2) is very close to the boundary of the T = 0 infrared 
instability region of the model (1.5). This behavior is rather generic and leads to the 
conclusion that the height size of the dome is very limited, the Tc is very small and 
not accessible through stable numerical analysis. Another way of realizing this issue 
relies on noticing that the BF bound is very mildly violated in the dome region such 
that the instability is very soft. 

We have repeated a similar dome analyses for the model (1.4), with the parameter 
A = am playing the role of disorder-strength in the boundary theory. Interestingly 
enough, we have found that again for A and q tuned to a small vicinity of the point 


25 












(7.2) we obtain a superconducting dome. This time, however, the normal phase can 
only be metallic. The superconducting dome at T = 0 is an interval [Ai, A 2 ], existing 
between two regions of normal metallic state, at A G [0, Ai] and A > A 2 . The critical 
temperature is bounded from above by a small number, and we did not access the 
finite-temperature superconducting state. We plot our results for the dome in model 
(1.4) in figure 14. 

This analysis shows that the existence of a superconducting dome region is a 
rather generic feature of these models, independent of the choice of the potential. 
In the next section we discuss the possible ways to alleviate the problem of the 
flatness of the dome. This seems to require introduction of an extra elements to our 
holographic system. 

8 Discussion 

In this paper we considered a holographic superconductor with broken translational 
symmetry, continuing the research, initiated in [5, 6]. To break the translational 
symmetry we used the known technique [12], coupling our system to the sector of 
massless neutral scalar fields, depending linearly on the spatial coordinates. We 
studied the standard Lagrangian for these neutral scalars, as well as its non-linear 
generalization, proposed in [17]. 

We have constructed models, exhibiting the following non-trivial new features: 

1. The Holographic superconductor in the non-linear Lagrangian model has a rich 
phase diagram on the temperature-disorder strength plane. In particular the 
superconducting phase is separated from the normal pseudo-insulating phase 
and the normal metallic phase by the line of second order phase transition as 
shown in figure 10. 

2. In the same model the optical conductivity exhibits a non-trivial emerging 
structure, signaling a collective excitation of the charge carriers localized in 
the mid-frequency range. This has been observed in [17] in the normal phase 
of the same model, for temperatures, lower than a certain critical value. In 
this paper we have demonstrated that this structure persists in the supercon¬ 
ducting phase. Eventually it gets destroyed by the charge condensate. This 
suggests a possible competition between the superconducting mechanism and 
the momentum dissipating one. In particular it seems clear that a large su¬ 
perfluid density completely screens this collective excitation which in a sense 
gets eaten by the large condensate. We are not aware of real superconducting 
system supporting a collective localized excitation like the one we see. In [17] 
this excitation was compared to a polaron excitation; it would be dehnitely 
interesting to make a comparison between the behaviour of this collective ex¬ 
citation in the holographic model and what really happens to a polaron when 
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superconductivity onsets. Unfortunately we are not aware of such a mecha¬ 
nism in real condensed matter systems. It would be also nice to see if other 
holographic models, providing translational symmetry breaking, support the 
same property. In this direction it would be very interesting to study the QNM 
structure of the system as initiated in [29]. 

3. We performed a complete analysis of the behavior of the critical (supercon¬ 
ducting) temperature as a function of the various parameters of our model. In 
particular we studied the curious non-monotonic behavour of Tc as a function of 
the graviton mass m, which was already observed in [5, 6]. Our results suggest 
that this feature persists for generic Lagrangian for neutral scalars. We do not 
have any clear explanation of the big mass regime where actually increases 
with the strength of translational symmetry breaking. It is even tempting to 
doubt the model in that regime, reminiscing the following known issues: for 
large momentum dissipation it seems that the energy density of the dual held 
theory at zero charge density gets negative [29]; the diffusion bounds for the 
model are unrestricted from below and the diffusion constants go to zero in 
that limit [30]. 

A very similar behavior has been observed in holographic SC with helical lat¬ 
tices [7] and with disorder [3]. It would be interesting to further analyze the 
universality and the meaning of this feature. 

4. By tuning the values of the scaling dimension A and the charge q of the scalar 
held, which is a bulk dual to the charge condensate of the boundary supercon¬ 
ductor, one can obtain a system, which exists in a superconducting phase, en¬ 
closed in a dome region. The dome region occurs upon increasing the disorder- 
strength parameter of the model, behaviour which is dehnitely diherent from 
the actual High-Tc SC phase diagram where this happens because of the dop¬ 
ing of the material. The critical temperature of the dome is very small, and 
in fact appears to be too hard to calculate numerically. The superconducting 
dome exists for both linear and non-linear models. In the case of the model 
with the linear Lagrangian for the neutral scalars, the superconducting dome 
exists in the middle of the normal metallic phase. In its non-linear extension 
instead, the dome exists between a pseudo-insulating phase for smaller values 
of disorder strength, and a metallic phase for larger values. There are no ex¬ 
perimental evidences of dome regions occurring because of disorder. We hope 
our work could in a way motivate some experimental effort in that direction. 

We are aware of two only holographic examples which show a superconducting 
dome region in way different setups [31, 32]. 

It would be interesting to improve the model, so that pseudo-insulating phase is 

replaced by an actual insulating phase. This will make the phase diagram more 
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resembling such of an actual high-Tc superconductor. This could be easily achieved 
introducing a dilaton held into the model. 

One direct expectation of our Holographic model is that the dome of supercon¬ 
ductivity we hnd seems to only exist in a very hne tuned region of parameter space 
which is always very close to the zero temperature instability. As a conseguence two 
immediate questions arise: 

• Is it possible to enlarge signihcantly the region of the parameter space where 
th dome appears? 

• Is it possible to get a dome with a reasonable which can be numerically be 
resolved? 

Solving the second issue would be indeed very important to rule out possible non-IR 
instabilities that would remove the dome region. 

One of the ways to accomplish this might be realized by the inclusion of a non¬ 
trivial coupling K between the charged scalar condensate and the neutral scalars as 
already shown in the action 2.2. Generically it seems that without the introduction 
of additional elements this can not be obtained. We leave this question for future 
work. 

A further interesting question is to look at universal properties of these large 
class of effective toy models such as the accomplishment of Homes’ Law following 
[33]. We leave these interesting questions for future investigation. 
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A Condensate and grand potential 

The aim of this appendix is to provide more details about the computations and the 
numerical procedures we did in Section 5. 



Figure 15. Condensate for A = 2, q = 1, auh = 0.5, rriL = 1 model with V = z + z^, 
and the corresponding imposed near-horizon data (for xi'^h ~ c) = !)• 

A.l Condensate 

In this subsection we will outline the routine to obtain the numerical solution of the 
equations of motion (2.7)-(2.10) for the whole superconducting background. First of 
all, evaluating the equations (2.7)-(2.10) at m = Uh, we can express 
x'{uh), A'liuh) in terms of 'ip{uh), xi'^h), A[{uh). Therefore we impose the initial 
conditions at Uh — e in the following way: 

-e) = ii{uh) - e'<p'{uh) , tp'iuh - e) = i)'{uh ), 
f{uh -e) = -ef{uh ), xiuh - e) = xM - ex\uh ), (A.l) 

e 2 

A{uh - e) = -eA[(uh) - —A”{uh) , A[{uh - e) = A^{uh) - eA''(uh) . 

where e is a small IR cutoff. One can solve equations of motion near the horizon 
to arbitrary order in e. We have found that imposing (A.l) is sufficient. We have 
checked explicitly that the results are stable towards changing e. We have the freedom 
of choice of the initial conditions 'ip{uh), A[{uh), and xAh)- The freedom of choice 
of xAh) is spurious, due to the time scaling symmetry, as we discuss below. 

The values of 'ip{uh) and A[{uh) are fixed by the requirement of having a fixed 
temperature T /and zero source 'ipi = 0, see (4.10). Both the charge density p, 
in units of which me measure the temperature, and the source ipi are determined by 
the near-boundary behavior of the numerical solution, with the gauge field behaving 
as: 

At{u) = p — pu + 0{A). (A. 2 ) 

In practical calculation we do the following. Suppose the temperature is suf- 
hciently small, so that the system is in a superconducting phase. We know that 
increasing the temperature will decrease the condensate, " 02 /^^^, until hnally at the 
critical temperature the condensate is zero. At that point 'ip{uh) = 0, that is, we 
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do not have the solution with vanishing source and non-trivial prohle of 'ip{u) in the 
bulk. Therefore we can start at ‘ip{uh) = 0, and take gradually incrementing values 
of 'ipiuh)- For each value of 'ip{uh) we search for A[{uh), such that 'ipi = 0. For an 
example of this kind of result see hgure 15. Finally for the given pair {'ip{uh), A[{uh)) 
we calculate numerically as shown for example in hgure 9. 

Scaling symmetry 

The equations of motion (2.7)-(2. 10) are invariant under the scaling symmetry: 

u = u/a, (t,x,y) = (t,x,y)/a, At = Ata, a = aa, (A.3) 

where a is a parameter of the symmetry transformation. The temperature, chemical 
potential, and the charge density therefore transform as: 

T = aT, yL = afl, p = p. (A.4) 

The scaling symmetry (A.3) allows one to hx Uh = 1. If Uh is not hxed to one, then 

we should substitute uf^‘^A[{uh) as the initial condition for the hux at the horizon. 

We have checked explicitly that the results are invariant under change of Uh- 

Time scaling symmetry 

The equations of motion (2.7)-(2.10) are invariant under the time scaling sym¬ 
metry: 

e^ = 6 ^e^, t = bi, At = At/b, (A.5) 

where 6 is a parameter of the symmetry transformation. 

We can use the time scaling symmetry (A.5) to £x x(0) = 0 at the boundary. 
This is necessary, so that the speed of light in the boundary held theory is equal to 
one. To achieve this, we impose the initial conditions on x to be xi.'^h) + 2 log b, and 
on the hux to be A[{uh)/b. We hx xi^h) once and for all. We have demonstrated 
explicitly that the result is independent of the choice of xi^h)- 

After hxing xi^h), for the given A^{uh), we integrate numerically the equations 
of motion, with 6 = 1. We then impose 6 = where x(0) is determined 

numerically. For this 6 we impose the initial conditions xi'^h) + 2 log 6 , A[{uh)/b and 
integrate the equations of motion again. This time, due to the time scaling symmetry 
(A.5), we have x(0) = 0. We have verified this explicitly. 

Running the described numerical procedure we were able to construct the con¬ 
densate ' 02 /p^A as a function of temperature hgure 9 we provide the plot 

of the condensate, for the model (1.5) with A = 2, q = 1, auh = 0.5 (a in units of 
entropy density), mL = 1. For the same parameters we also plot the initial condi¬ 
tions {—A'^{uh),'ip{uh)), which we imposed, to enable the vanishing source = 0 in 
hgure 15. 
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A.2 Grand potential 

Here we provide intermediate steps for calculation of the grand potential. 

The holographic prescription for the calculation of the grand potential is: 

Vt = -T \ogZ = TSe, (A.6) 

where is a Euclidean on-shell action of the bulk theory. This should be supple¬ 
mented with the boundary Gibbons-Hawking term, and the counter-terms^®. The 
resulting action reads: 

Se = I + Igh + Ic.t ., (A.7) 

where the boundary Gibbons-Hawking term is given by: 


Igh = -2 


d^x \/—hK 




U=€ 


(A.8) 


where e is a UV cutoff, hab the pullback metric on the boundary and Kah the extrinsic 
curvature^^. The counter-term action I^.t. is a sum of gravitational, scalar and axion 
helds counter-terms [5, 34]: 


Lf = 


d^x 


Y + y — 2m^ L V 

Jj 1 j 


It is convenient to evaluate the following Lagrangian on shell: 

L = L + 2 (y^K^ . 

to get: 

I + Igh = j d^x L — 2 J d^x y/—h K 

U=Uh 

After a straightforward calculation we obtain: 


(A. 10) 


(A.11) 


(A. 12) 


(A. 13) 


Notice that which B{uh) = 0. Therefore the full on-shell action is given by: 

I + Igh = j d^x - 5(e)) • (A.14) 

^^Equivalently, one can calculate difference of the grand potentials of two phases, in order to 
avoid adding the counter-terms. 

^^It is defined by 

K = V^n^, = (O, 0, 0, (A.9) 

where is the unit vector normal to the boundary. 
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To proceed with the calculation, we need to be able to evaluate the counter-term 
action (A. 10) and the B{e) term of (A. 14). We need to know the near-boundary 
behavior of the fields. That is given by : 


At = ^ — pu + 0{u‘^) 


(A.15) 



(A.16) 

(A.17) 

(A.18) 

(A.19) 


/ — 1 + 7iU^ + 72U^ + 0(u‘^) 

X = + Gv? + 0{u‘) 

V = Viu‘ + 0(11®). 


We are interested in the systems with vanishing source of the charged scalar, 'ipi = 0. 
By solving equations of motion near the boundary, we obtain: 


7 i = Vi. 


(A.20) 


Combining all the results together, we arrive at the hnal expression for the on-shell 
action; 



(A.21) 


B On-shell action for fluctuations 

The calculation of the on-shell action for fluctuations is similar to the one for the 
grand potential performed in appendix A. The total action is a sum of the total 
bulk action (2.1), the Gibbon-Hawking (GH) term on the boundary (A. 8), and the 
counter-term action (A. 10): 



(B.l) 


We evaluate the action (B.l) on the ansatz: 



^—f{u)dt‘^ + 2ehtx{u, t)dtdx + dx^ + dy^ + 



A = Afdt -|- e ax{u, t)dx , 
(ff = ax + e^iuA), 

(l)y = ay, 

V’ = '0(m) , 


(B.2) 


and collect C>(e^) terms, which describe dynamics of the fluctuations htx-, o-x, C- The 
0(e) terms vanish due to equations of motion, satisfied by the background fields /, 

^®Note that this is true only if the potential reads V (X) = X + -|-... where the smallest 

power is always equal to one. 
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At, (j)^'y, 'ip, and the 0{e^) terms are contributions to the grand potential for the 
background. 

The GH term vanishes at the horizon. Therefore: 

!’ = ll + Igh = h - (-2^A')' , (B.3) 

We obtain: 

(-2/e^ {2uhtx{u,t) {u^fA[duax{u,t) + 2q^uAt'ip‘^ax{u,t) 

+ Afduhtxiu, t) + 2aL‘^w?udti{u, t)v'j + (- {u^{dtax{u, t)Y 

+ f{dthtx{u,t)f + 2L^m^{dti{u,t)fvy^ + htx{u,tf {-2{uf' + 2,) 

+ / + 2nx' - 6) + {2m^ {V - a^u^V) + 

- u^e^^htx{u,tf {u^A't + 2q^Al'ip‘^) - 2u^ f{duax{u,t)f 

- Aq^u'^ftp'^axiu, - 4L^mV/^(5„^(n, t)) Vj (B.4) 


To proceed, we integrate the a'^, terms by parts, and substitute expressions 

for a", h'!^, from the corresponding fluctuation equations. We need to keep track 
of the boundary terms. Then let us go to the momentum space. As a result we arrive 
at If = B'j, where: 


Bf = 


L^e-t 

2m3 


(e^htxiu, -uj) {u^A[ax{u, u) - uh't^{u, oj) + Ahtx{u, cn)) 


+ uf (u^ax(u, —oj)a'^{u, co) + 2L?m^^{u, —oj)lP{u, oj)V 


(B.5) 


where prime, as before, stands for a derivative w.r.t. u. 

The counter-term action (A.10) for the ansatz (B.2) is given by: 

_ X 

lit = 2u^p/‘i {htx{u, -u:)htx{u,u:) (2LW1/ 

— 'ip‘^—4:) —2L‘^m'^u^V (^u‘^^{u,—oj)^{u,oj)+ahtx{u,—oj){ahtx{u,u))+2iu)^{u,u)))'^'^ , 

(B.6) 


evaluated at u = 0. 

Now let us evaluate (B.6) minus (B.5) at n = 0, which gives ll^t- Consider the 
case A = 2. First we need to solve fluctuation equations near the boundary. We 
already determined the near-boundary asymptotics (A. 17) for the background helds. 
In superconducting phase we have 'ipi = Q. Besides, from the equations of motion, 
one obtains: 


T / ‘1 T ‘1 2 

7 i = —Vi m L a . 


(B.7) 





Similarly, the fluctuation equations of motion, near the boundary give: 


ax{u) = u , (B.8) 

htx{u) = h[l\uj) + h[l\uj)u^ + hll\uj) , (B.9) 


where again not all the coefficients of expansion are independent, and in fact: 

= 2Vim‘^apu 

hil\uj) = Yp (2^1 + o,x\^) (B.12) 

hfx{^) = Y “ 6iVim^a^^^\u)) . (B.13) 

Using these asymptotic expansions, evaluating l[^ — Bf ai u = 0 gives^® : 

iL = ax\-(^)ax\^) - pax^hil\-uj) + 272 h[l\-uj)hil\u) - 3hll\-u)hil\u) 

+ 6m‘^Vi^^^\-u)^^^\u). (B.14) 


where we have kept for brevity (but keep in mind it is not an independent 
expansion coefficient, dne to (B.IO)). 

It is convenient to replace ^ Z, so that we are dealing with two fields, {ax , Z), 
which have the same near-bonndary expansion, at least np to the first two orders. 
Dne to (6.5), we obtain: 

Z{u,u) = , (B.15) 

lojau 

which near the boundary becomes: 

Z{u, oj) = Z^~^\oj) + Z'^‘^\oj) « + ••■ =-- -- ^*'^^(n;) n + ... . (B.16) 

auj auj 

We can represent in the form, convenient for calculation of correlation matrix: 

IL = (4‘'(-‘^). M j + ■ ■ ■ . (B.17) 


where dots denote terms. We cannot extract by solving system of eqnations 
for {ax, Z), because Z ~ So we assume that is a constant of integration, 
which we fix to be: 




i(l + \/2){u}‘^ — 2m?a‘^Vi) 
2m?apuVi 



(B.18) 


^®This is in agreement with eq. (3.14) of [6]. See that only the leading linear term Vi in the 
near-boundary expansion of V (z) matters in this formula. 
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which is the choice enabling a diagonal matrix M. Let us rescale the fluctuation 
flelds (this is a symmetry transformation of fluctuation equations): 


-)■ 


The corresponding matrix is: 


(B.19) 


M = 


( \ 

1 0 

g __ 

^l-2V2+^^^ ) 


(B.20) 


For the purpose of finding AC conductivity we only need the a^) component of 
the correlation matrix. 
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